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Abstract 

It is shown that the so-called “Yukawaon” model can give a unified description of masses, 
mixing and CP violation parameters of quarks and leptons without using any hierarchical 
(family number-dependent) parameters besides the charged lepton masses. Here, we have 
introduced a phase matrix P = daig(e l<)il , e 1 ^ 2 , e 1 ^ 3 ) with the phase parameters (fi, <j>i, ff) 
which are described in terms of family number-independent parameters, together with using 
only the charged lepton mass parameters as the family number-dependent parameters. In this 
paper, the CP violating phase parameters Sf, p and Y CP in the standard expression of Vckm 
and Upmns are predicted as 8 q CP ~ 72° and 6 CP ~ — 76°, respectively, i.e. S e C p ~ ~5 q c P - 

PCAC numbers: 11.30.Hv, 12.15.Ff, 14.60.Pq, 12.60.-i, 

1 Introduction 

1.1 What is the Yukawaon model 

One of the big subjects in the particle physics is to investigate the origin of flavors. There 
is an attractive idea that the flavor physics is understood from the point of view of a family 
symmetry [T]. However, the symmetry has to be explicitly broken by the Yukawa coupling 
constants Yj- (/ = u, e, u , d) if we suppose the family symmetry to be a continuous symmetry. 
Therefore, the symmetry is usually considered as a discrete symmetry. If we adhere to the basic 
idea that the flavor symmetry should be a continuous symmetry which is unbroken at the start, 
we are forced to consider that the Yukawa coupling constants are effective coupling constants 
Yj^ which are given by vacuum expectation values (VEVs) of scalars (“Yukawaons”) Yj with 
3x3 components [2]: 

{ Y f ff )i 3 = (/ = u,d,v,e ), (1.1) 

where A is an energy scale of the effective theory. In the Yukawaon model, all the flavons [3] are 
expressed by 3 x 3 components of U(3). We consider no substructures of it such as 2 x 2 and so 
on. 
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In the Yukawaon model, we assume a U(3) family symmetry, and take the following would- 
be Yukawa interactions: 

Hy = + ^(hnYeKMHd + yR{p R ) i {Y R ) ij {v c R y 

+^(q L ) i (Y u ) i j (u R ) j H u + ^(qLYiY^idnyHd, (1.2) 

where £l = (zy+l) and q R = ( ul,cIl ) are SU(2 )r doublets. H u and Hd are two Higgs 
doublets. The third term in Eq.(1.2) leads to the so-called neutrino seesaw mass matrix [@] 
M u = YyY^Yj , where Y v and Yr correspond to the Dirac and Majorana mass matrices of 
neutrinos, respectively. Hereafter, for convenience, we use notation A, A and A for fields with 
8 + 1,6 and 6* of U(3), respectively. 

In order to distinguish each Yukawaon from the others, we assume that Yf have different R 
charges from each other by considering P-charge conservation [a global U(l) symmetry in N = 1 
supersymmetry (SUSY)]. Of course, the P-charge conservation is broken at an energy scale A, 
at which the U(3) family symmetry is broken. For R parity assignments, we inherit those in the 
standard SUSY model, i.e. R parities of yukawaons Yf (and all flavons) are the same as those 
of Higgs particles (i.e. Pr( fermion) = —1 and PR(scalar) = +1), while quarks and leptons are 
assigned to Pr (fermion) = +1 and Pr (scalar) = — 1. 

A remarkable characteristic of the Yukawaon model is that it is possible to understand 
the observed hierarchical structures of masses and mixings of quarks and leptons without using 
any family number-dependent parameters except for the charged lepton masses. That is, all 
the quark and lepton masses and mixings can be understood in terms of only the observed 
charged lepton masses. In the Yukawaon model so far, our aim seems to be almost accomplished 
except for the following problem only: Namely, we have obliged to introduce a phase matrix 
(P) = diag(e*' ?il , e*^ 2 , e*^ 3 ) in order to give a good fitting for the Cabibbo-Kobayashi-Maskawa 
(CKM) mixing matrix [5], where (fa, fa, fa) have been introduced as family number-dependent 
parameters. In this paper, however, we will relate those parameters (4> 1, fa, fa) to the observed 
charged lepton masses m e i as discussed in Sec.4. 

Relations among Yukawaon VEVs (Yf) are obtained by supersymmetric vacuum conditions 
from U(3) symmetric and P-charge conserved superpotential. In this article, the Yukawaon 
VEVs (Yf) are related to VEVs of fundamental flavons <f+ with a common bilinear form to all 
flavors o 

(Y f ) = k f (<S> f )(H> f )+fal C f = u,d,u,e ), (1.3) 

1 In earlier Yukawaon models [ 6 ], the bilinear form was only for the up-quark sector, and the model could 

excellently lead to the so-called tribimaximal mixing [7] in the lepton mixing under the use of only a few param¬ 

eters. However, the model [ 8 ] (and also [ 8 ]) could not give the observed sizable mixing angle #13. We found that 
we can give the observed value sin 2 2+3 ~ 0.1 [5] only when we consider that all VEV relations are given by a 
common bilinear form Eq.(1.3) [10] . 
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where 1 = diag(l,l,l). Here, the VEV matrices (Ty) and (Tj) are commonly related to a 
fundamental flavon VEVs (To) and (‘ho) by 


(Pf) ik (*f)ki(P f ) lj = k' f (^(s f U(^, 

(Pf^f^iPfh = k'f(®o)i a (S f ) a P ($o) /3j, 


(1.4) 


and Pf 


(1.5) 


( 1 . 6 ) 

The form of Eq.(1.6) is understood by a symmetry breaking U(3) / —>- 83 . 

1.2 Charged lepton sector as a fundamental flavor basis 

We consider that the charged lepton mass matrix is the most fundamental one compared 
with other mass matrices and that the charged lepton mass values play an essential role in 
understanding the flavor physics. The points of our postulation are as follows: 

(i) There is a fundamental flavon To, and the reference basis in the flavor physics is defined by 
the diagonal basis of (To) an d (To): 

(T 0 ) = (T 0 ) = u 0 diag(xi,x 2 ,x 3 ). (1.7) 

where x* are real parameters with x\ + x\ + X 3 = 1 . 

(ii) In the reference basis, the U(3) / family symmetry is broken into S 3 , i.e. VEVs of flavons Sf 
and Sf take the form (1.5). 

(iii) The charged lepton mass matrix ( Y e ) should be diagonal and real as well as (To) and (To), 
and it should be described only in terms of the fundamental parameters Xj. Therefore, with 
demanding simplicity too, we require 

a e = 0, £e = 0. (1.8) 

This means Xi oc rri l J 4 ( m e i = (m e , m T )). In Sec.3, we use the following parameter values of 

Xi, 

(xi,x 2 ,x 3 ) = (0.115144,0.438873,0.891141). (1.9) 

2 The form (1.5) was suggested by a “democratic universal seesaw” mass matrix model [11], in which quark 
mass matrices are given by a form (4> e )(l + a/.Y3)(4> e ). 


where i and a are indices of U(3) and U(3)', respectively. For the VEV structures of Pf 
in Eq.(1.4) and the £/ terms in Eq.(1.3), we discuss in the next section. 

The VEV structures (Sf) and (Sf) in Eq.(1.4) are given by! 


where 


(Sf) = (Sf) = 1 + afX 3 , 
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In (1.9), we have used running mass values m e ([x) = 0.000486847 GeV, m^/i) = 0.102751 
GeV and m T (n ) = 1.7467 GeV as the charged lepton mass values at /j = Mz , because our 
numerical predictions in the quark mass ratios are done at p, = Mz- Note that the mass values 
(? 7 i e , m r ) have a large hierarchical structure, i.e. m e /m T ~ ICG 4 , while the values (1.9) have 
a mild hierarchical structure, i.e. x\/x 3 ~ 10 ~ . 

In this paper, we do not ask the origin of the value (xi,X 2 ,xs). However, for reference, in 
Appendices A and B, we will demonstrate an example of the charged lepton mass relation in 
the present Yukawaon model. 

1.3 What is new 

The new characteristic points of the present Yukawaon model are as follows: 

(i) The bilinear form of the Yukawaon VEVs given in Eqs.(1.3) and (1.4) has been adopted for 
some flavor sectors in the previous paper uni, too. However, in the present paper, we apply it 
to all the sectors f = u,d, v, e. 

(ii) So far, we have needed a VEV matrix P given by 

P = diag(e^ 1 ,e^ 2 ,e^ 3 ), ( 1 . 10 ) 

in order to fit CKM mixing matrix Vckm reasonably. However, our aim was to describe all the 
masses and mixing of quarks and leptons in terms of family number-independent parameters 
expect for the charged lepton masses. Therefore, the phase parameters (f>i were against our aim 
and unwelcome as it is. In the present paper, we try to denote the family number-dependent 
parameters cf>i in terms of the observed charged lepton masses rn cl . The details are discussed in 
Sec.4. 

(iii) In general, the phase matrix P affects not only the CKM quark mixing matrix Vckm 
but also the Pontecorvo-Maki-Nakagawa-Sakata [T2] (PMNS) lepton mixing matrix Upmns- 
Furthermore, predicted value of the CP violating phase parameter 8 £ CP in Upmns will depend 
on an appearing position of P which is arbitrary at present. In the present paper, as we will 
discuss in Sec.2, we construct a model in which the phase matrix P appears in the up-quark 
sector only, and it affects not only the Vckm but also the Upmns through Eq.(2.5) as we give 
later. As a result, as we discuss in Sec.3, we predict CP violation phase parameters 5 q cp and 5 l CP 
in the standard expression of Vckm and Upmns as 5^ P — —8 q CP — —70° unlike the previous 
papers. 

In Sec.2, we construct a mass matrix model base on a Yukawaon model in which the phase 
matrix P appears in the up-quark sector. In Sec.3, parameter fitting are discussed. Especially, it 
is shown that CP violating phase parameters 5 q CP and 8q P in the standard expression of Vckm 
and Upmns are predicted as 5 q CP — 72° and 8 £ CP ~ —76°, respectively, i.e. Sq P ~ ~^cp- I n 
Sec.4, we will propose a new relation between P and m e i- Finally, Sec.5 is devoted to concluding 
remarks. In Appendices A and B, we will demonstrate an example of the charged lepton mass 
relation in the present Yukawaon model. 


4 


2 Model 


2.1 VEV relations and superpotentials 

The VEV relations of Yukawaons have already been given in Eqs.(1.3) and (1.4) in Sec.l. 
Those VEV relations are derived from a U(3)xU(3) / symmetric superpotential. For example, 
for the VEV relation (1.3), we assume [10] 


Wy= V {tt [(jifYf + v®/®/) ®f] + Tr [m/P + *'/*/*/] T 1 '!®/]} • (21) 


f=u,d,v,e 


The supersymmetric vacuum condition dW/dOf = 0 leads to 


Vf Y f + ^ [^f Y f + 

that is, to the VEV relation (1.3) with the £/ term, 


1 = 0 , 


( 2 . 2 ) 


£/ 


^/ A / ~ Vf X 'f 
(Pf + 3 


Tr[d>/4>/]. 


(2.3) 


Here, we assume that flavons Of always take (0/) = 0, so that vacuum conditions, which are 
obtained by differentiating the superpotential IVy with respect to other flavons, do not give any 
additional VEV relations, because those relations always include one (0/). 

In a similar way, we can obtain the VEV relation (1.4) by assuming 

^ = E ( A /1 (Pf) ik (*fMPf) lj + *f2($>o) ia (SfU($o)P J ) (0/)ii, (2.4) 

f=u,d,u,e 


and so on. Although we assumed a tad pole term Tr[0j] in (2.1), we consider that such a term 
is a special case only for Of with (8 + 1) of U(3). Therefore, such terms do not appear in 
the relation (1.4). 

In order to give a neutrino mass matrix with seesaw mechanism, M u = k l ,(Yj')(Yf>)~ 1 (Y u ), 
correspondingly to the following Majorana neutrino mass matrix M }j (i.e. (Yr)), we assume 
VEV matrix relation 

(Y R ) = (Y e )(<S> u ) + (<Z> u )(Y e T ), (2.5) 

according to the previous Yukawaon model m- 

2.2 VEV structures of Pf 

Prior to discussing VEV forms (Pf) and (Pf) given in Eq.(1.4), let us consider the following 
superpotential 

Wp = Etc [PPEE] + ^-Tr[PP]Tr[EE], (2.6) 
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where, in order to distinguish P from E, we assign R charges of P and E as 


R{P) = R(P) = ^(1 - A), R(E) = R(E) = i(l + A), (2.7) 

so that R,(P) + R(P) + R{E) + R(E) = 2. The supersymmetric vacuum conditions lead to 

(P)(P) = 1, (E)(E) = 1 . (2.8) 

We dehne specihc solutions of (2.8) as 

(P) = diag(e^,e^,e^), (E) = diag(l,l,l). (2.9) 

We consider that VEV of each flavon ( Pf ) given in Eq.(1.4) is given by either ( P ) or (E) in 
Eq.(2.9) under the 77-term condition as discussed in (2.10) below. 

On the other hand, let us recall that, in general, VEV matrix (yl) is related to VEV matrix 
(A) under the D term condition as 

(A) = (A)*, or (A) = (A). ( 2 . 10 ) 

Let us back to the relations (1.4). We take 

(4* 0 ) = (4>o) = diag(xi,x 2 .x 3 ), (2.11) 

(S f ) = (S f ) = l + a f e ia fX 3 , (2.12) 

while 

(. P ) = (P)* = diag(e-^,e-^,e-^), (2.13) 

where parameters xy, fa, a,f and af are real. Then, according as (4> f) = (<3?j)* or (4>j) = ($f), 
the relations (1.4) require (Pf) = (Pf)* or (Pf) = (Pf). For example, when we take the case 
(!>/■) = (4*/)*, Eq.(1.4) becomes 


(Pf)*(<S> f )(P f )* = (* 0 )(Sf)(* o), 
(P f )($ f )*(P f ) = ($o)(S f )($ 0 ). 


Two equations in (2.14) cannot simultaneously satisfied without ctf = 0. On the other hand, 
when (4>j) = (4>j), Eq.(1.4) becomes 


(PfY(*f)(PfY = ($o)(5/)(4>o), 
(Pf)($f)(Pf) = ($ Q )(S f )($o). 
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Then, for the case of aj fa 0, the two equations (2.15) are satisfied only when fa = 0 [i.e. 
(Pf) = {E)]. As a result, we consider only two cases: (i) for the case of aj fa 0, we regard (Pf) 
as (E) given in (2.9), and (ii) only for the case of atf = 0, we regard (Pf) as (P) with fa fa 0. 

The parameters aj affect not only CP violation, but also mass ratios. In the up-quark 
sector, as we discuss in Sec.3, we can fit up-quark mass ratios m u /m c and m c /mt by taking two 
parameters a u and fa (keeping a u = 0). Therefore, we regard the up-quark sector as the case 
of OLf = 0, so that we regard (Pf) as (Pf) = (P). Also, since we have taken a e = 0, we have to 
regard (P e ) as (P e ) = (P). Note that (P) and (Y e ) are diagonal, so that they are commutable 
each other. Therefore, (P e ) does not play any essential physical role in the parameter fitting of 
the masses and mixing of quarks. Hereafter, we denote (P e ) as (E) from the practical point of 
view, except for a case of counting of R charge. On the other hand, in down-quark sector, we 
cannot fit down-quark mass ratios rrid/m s and m s /mi) without help of ad fa 0. Therefore, we 
regard down-quark sector as a case of (Pf) = (E). Thus, we have the selection rule, (Pf) = (P) 
or (Pf) = (E), as a phenomenological one. For neutrino sector, we have no phenomenological 
information. For simplicity, we take a fewer parameter scheme (a u fa 0 rather than (fa fa 0). 
Hereafter, we will use the notation Pf as 


Pf = P for / = u, e, 
Pf = E for / = d, v. 


(2.16) 


Sometimes, for convenience, we use notations P u , P e , and so on, although we identify P u and 
P e as one flavon P , and also Pd and P u as one flavon E. 

The phase matrix (P) does not affect mass ratios. Since (P u ) = (P), the phase parameters 
affects CP violation phase 5)Lp in the CKM mixing matrix Vqkm ■ However, note that the 
phase in (P u ) = ( P ) can also affect CP violation phase 5q P in the PMNS mixing matrix 
Upm ns t because the phase in P u can affect Yp through as shown in Eq.(2.5). This is the 
most notable point in the present paper§| 

The details are discussed in the next section, Sec.3. 

2.3 R charge assignments 

In this model, the number of flavons is larger than that of VEV relations. Therefore, in 
general, we cannot uniquely determine R charges of flavons. Since we demand to assign R 
charges as simple as possible, we put the following rules for simplicity: 

(i) We assign the same R charge to flavons A and A: 

R(A) = R(A), (2.17) 

3 A similar model with P u has been discussed in Ref.JTO]. However, we has neglected a possible effect of { P u 
in the neutrino mixing. 
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independently whether (A) = (A)* or (A) = (A). Then, we obtain R charge relations 


R(Y f ) = 2R(^ f ) = 2r f (f = u,d, u,e), (2.18) 

and 

R($ f ) = R($ f )=R(S f ) + 2R($ 0 )-2R(P f ) (/ = u,d,v,e), (2.19) 

from Eqs.(1.3) and (1.4), and 


R(Pu) = R(Pe) = R{P) = \{l + A), 

R(Pd) = RiPu) = R{E) = i(l - A), 

from Eqs.(2.16) and (2.7). Therefore, from Eq.(2.19), we obtain the following relations: 


R(*e) - R(S e ) = R($ u ) - R(S U ) = 2R(4> 0 ) - (1 + A), 
R(* v ) ~ R(S V ) = R($ d ) - R(S d ) = 2i?($ 0 ) - (1 - A). 


(ii) We can regard that R charges of Yj are determined only by those of the SU(2)^ singlet 
fermions f c . Therefore, we simply assign 

R(£H U ) = R{£H d ) = R(qH u ) = R(qH d ) = r H + 2. (2.22) 

Since those have different quantum number of U(l)y, we can distinguish those from each other 
in spite of the relation (2.22). Then, we obtain a simple R charge relation 

R(Y f ) + R(f c ) = -r H . (2.23) 

For Yfi, we obtain 

R(Y r ) = 2 - 2 R{u c ) = 2 r H + 2 + 2 R(Y V ), (2.24) 

from Eqs.(2.1) and (2.23). On the other hand, from Eq.(2.5), R{Yr) must be satisfied a relation 

R(Y R ) = R($ u ) + 2R(<i> e ). (2.25) 

From Eqs.(2.24) and (2.25), we have the following constraint 

2 R($ e ) - 4 R(§„) + R{$ u ) = 2 r H + 2. (2.26) 

Even under the these constraints, we cannot still completely fix the R charges of whole 
flavons. In the present model, R charge assignments are not so essential, so that it is enough 
to assign R charges to distinguish flavons with the same U(3) from each other. That is, we are 
satisfied with any i?-charge numbers which satisfy the relations (2.18) - (2.26). Nevertheless, it 



is desirable to have explicit i?-charge assignments as simple as possible. Therefore, let us go on 
our search for explicit i?-charge assignments. 

First, for simplicity, we put 

m o) = (2-27) 

Then. Eq.(2.21) becomes to be simpler relations 

R(Sf) = R($ f ) + A (f = e,n), 

R(S f ) = R($ f ) — A (f = u,d). 

Now, let us discuss possible i?-charge assignments for Yukawaons Yf under the conditions 
discussed above. If we have R{Yf ) = 0 , then we can attach the field Yf on any term in su¬ 
perpotential. Therefore, we require R{Yf) / 0 for any f = e,v, d, u. Also, we have to require 
RiYfYf) ^ 0 for any combination of / and f. As a result, we have to consider that whole 
R values of Yf are positive. Furthermore, we speculate that the values of R will be describe 
by simple integers. Of course, the R charges have to satisfy the relation (2.26). Therefore, we 
assign simpler R charges to the Yukawaons Yf on trial as follows: 

( R(Y e ),R(Yu),R(Y u ), R(Y d )) = (1, 2,3,4), (2.29) 

that is, 

(R($ e ),R($ u ), R($„), R($ d )) = Q, 1, |,2 ) . (2.30) 

This assignment satisfies the condition (2.25) for R(Yr) with rn = —3. 

In Table 1 , as a summary of Sec.2, we present the assignments of SU(2)x,xSU(3) c xU(3)xU(3) / 
and the R charges of the fields in the present model. 

3 Parameter fitting 

3.1 How many parameters? 

Our mass matrices Yf for / = e, v. d, u with the VEV relations discussed in Sec.l and Sec.2 
are summarized as follows: 


Y e = $ e $:, 

= P*$ 0 <I> 0 P*, 


= PTqTqP 
d’o = diag(xi,x 2 ,x 3 ), 


(3.1) 
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Table 1: Assignments of SU(2 )lxSU(3) c xU(3)xU( 3)'. For R charges, see subsection 2.3. We 
assign the same R charges for flavons A and A, e.g. R(A ) = R(A). For a special choice, r e , r u , 
r u and r d are taken as r e = 1/2, r u = 3/2, r u = 2/2 and r d = 4/2. 



i = (is, e ) f c = is c , e c 

q = (■ u, d) f c = u c , d c 

H u H d 

SU(2) l 

SU(3) C 

2 1 

1 1 

2 1 

3 3* 

2 2 

1 1 

U(3) 

U(3)' 

3 3* 

1 1 

3 3* 

1 1 

1 1 

1 1 

R 

2 -(2r/ + r#) 

2 -(2 rf + m) 

r H ru 


% Yr 


+i 

O 

O 

^e,u ^e,u 

Su,d Sp,d 

1 1 

1 1 

1 1 

1 1 

1 1 

1 1 

i i 

i i 

1 1 

1 1 

8 + 1 6 

1 1 

6 6* 

1 1 

3 3* 

3 3* 

i i 

6 6* 

1 1 

6 6* 

2 r f r R 

r f 

1/2 

r e ,u ~ A 

r v ,d + A 


P 

P 

E 

E 

©/ 

©ij 

0 $/ 

0 $/ 

1 

1 

1 

1 

l 

1 

1 

1 

1 

1 

1 

1 

l 

1 

1 

1 

6 

6 * 

6 

6 * 

8 + 1 

6 * 

6 

6 * 

1 

1 

1 

1 

1 

1 

1 

1 

+ A) 

1(1 

-A) 

2-2 r/ 

2 - r R 

1 - 

R(S f ) 
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Y v — V^V T £l/l> 


(3.2) 


$V = E $o(l + a u e ia ''X 3 )^ 0 E, 

Y u = $«$* +£ u l, 

= P*$ 0 (l + a u X 3 )<S> 0 P*, 
= P$ 0 (l + a u X 3 )$ 0 P, 


Y d = + Cdl) 

§d = E $0(1 + a d e tad X 3 )$ 0 E, 
Neutrino mass matrix with seesaw mechanism is given by 

M„ = Y V Y~ X Y V , 

Yr = Y e $ u + $ u Y e , 


(3.3) 


(3.4) 


(3.5) 


Note that 4?* in Eq.(3.1) and 4>* in Eq.(3.3) are not <3? e and 4> n , respectively. Here, for conve¬ 
nience, we have dropped the notations “(” and We also make no distinction of property 
under U(3)xU(3) / , i.e. we denote A and also A as A simply. Since we are interested only 
in the mass ratios and mixings, we use dimensionless expressions $0 = diag(aq, X 2 , x 3 ) (with 
x\ + x\ + x\ = 1), P = diag(e* ?il , e*^* 2 ,1), and E = 1 = diag(l, 1,1). Therefore, the parameters 
a e , a u , a u , ad, and £ d are re-defined by Eqs.(3.1)-(3.4). 

In the phase matrix P defined by Eq.(l.lO), physical values are only differences among 
(4> 1 , 4> 2 , 4>z)i so that we can take one of 4>i (i = 1 ,2, 3) as zero in the parameter fitting for Vckm- 
In this paper, we put (j) 3 = 0, so that free parameters are ^ 2 ). Note that, as we stated in 
Sec.2.2, P and P* in Eq.(3.1) do not affect Y e practically, because 4>o and Y e are diagonal, so 
that P and P* are commutable with $0 and Y e . 

Therefore, in the present model shown in Eqs.(3.1) - (3.5), except for the parameters 
{xi,x 2 ,x 3 ), we have 10 adjustable parameters, (a u ,a u ,£ v ), {a u ,£ u ), (a d ,a d ,£ d ), and (^ 1 ,^ 2 ) 
for the 16 observable quantities (6 mass ratios in the up-quark, down-quark, and neutrino sec¬ 
tors, four CKM mixing parameters, and 4+2 PMNS mixing parameters). Especially, quark mass 
matrices M u = Y u and M d = Y d are fixed by two parameters (a u ,^ u ) and (a d ,a d ,Cd), respec¬ 
tively. Note that those parameters are family number-independent parameters. Therefore, in 
order to fix those parameters, we use two inputs values, up-quark mass ratios ( m u /m c ,m c /mt ) 
and down-quark mass ratios (m d /m s ,m s /mb), respectively, as we discuss in the next subsection 
3.2. After the parameters (a u ,^ u ) and (a d ,a d ,^ d ) have been fixed by the observed quark mass 
rations, we have five parameters (a u , a u , £„) and (<j>i, ^> 2 ) as remaining free parameters. Processes 
for fitting those five parameters are listed in Table 2. In subsection 3.3, we discuss the fitting of 
four CKM mixing parameters, |V^ S |, \V c b\, \V u b\ and \V td \, by adjusting two parameters (0i,<fe)- 
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Table 2: Process for fitting parameters. N parameter and Ni nput denote a number of free parame¬ 
ters in the model and a number of observed values which are used as inputs in order to fix these 
free parameters, respectively. E N - means E N parameter or E N input 


Step 

Inputs Ninput 

P ar ameters Nparameter 

Predictions 

1 st 

2 nd 

3rd 

option 

m u /m c , m c /m t 2 

m d /m s , m s /m b 2 

\Vd,\, \V ub \ 2 

sin 2 2012, sin 2 20 2 3 , Ru 3 

A m§ 2 

Q"ii, iu 2 

a d , cud, id 3 

( 01 ,02 ) 2 

+ i Olu, iu 3 

n+3 

\Vus\, \Vtd\, s q CP 
sin 2 2013 , Sq P 

2 Majorana phases, ^ 

(mui,mu2,mu 3 ), (m) 

EN.. 

9 

10 



Also, in subsection 3.4, we do the fitting of PMNS mixing (sin 2 20i2, sin 2 2023, and sin 2 2013) 
and neutrino mass ratio ( R u = by adjusting three parameters (a v ,« v ,^). 

3.2 Quark mass ratios 

First let us fix values of (a u ,i u ) from the up-quark mass ratios. The observed values of the 
up-quark masses at ^ = mz ESI are 


r 


u _ 
12 — 


m v 


m. 


= 0.045 


+0.013 

- 0 . 010 ’ 



0.060 ± 0.005. 


(3.6) 


We obtain four solutions of (a u ,i u ) which can give the values (3.6). Among them only one 
solution 

K,e u ) = (-1.4715,-0.001521), (3.7) 

can give a reasonable prediction of the PMNS mixing as we discuss later. 

Secondly, let us fix values of (a d , a d , id) from the down-quark mass ratios. From the observed 
down-quark mass ratios m 

rf 2 = ^ = 0.0531°;°°*, 4 = ^ = 0.019 ± 0.006, (3.8) 


or [H] 


d - m d 

r 12 = - 

m .o 


= 0.050 


+0.002 

- 0 . 001 ’ 


4 = El = 0.031 ± 0.005, 
m b 


(3.9) 


we determine the parameters (ad, id, oi d ) as follows: 

(a d , OL d , id) = (-1.4733,15.694°,+0.004015), 


(3.10) 
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which leads to the numerical results as follows: rf 2 = 0.0612, r 23 = 0.0312. These values are 
inconsistent with the observed values (3.9), but, roughly speaking, those are consistent with 
(3.10). We think that the light quark mass values are still controversial. 

Here, we have tried to fix the parameters (ad,ctdi£d) in the down-quark sector by using 
input parameters [13] for rf 2 and r 23 . However, since we have three parameters for two input 
values md/m s and m s /rrib , we cannot fix our three parameters. It is more embarrassing that 
there is no solution of m s /mb ~ 0.019 in the (ad,ad,£d) parameter region. Nevertheless, we 
found that the minimal value of m s /mb is m s /m(, ~ 0.03 at (a^a^, £d) ~ (—1.47,16°,0.004) 
which can give a reasonable value of md/m s at the same time too. Therefore, we take the values 
in Eq. (3.10), which leads to rf 2 = 0.0612 and r 23 = 0.0312. Note that the value r 23 = 0.0312 is 
considerably large compared with r 23 ~ 0.019 by Xing et al. [13], while the value is consistent 
with r 23 — 0.031 by Fusaoka and Koide [Hj. The values md{y) and m s (fi) are estimated at a 
lower energy scale, fj, ~ 1 GeV, so that we consider that the ratio rf 2 at fi = Mz is reliable. On 
the other hand, the value is extracted at a different energy scale /x ~ 4 GeV from /x ~ 1 

GeV, so that the value rrib(Mz ) is affected by the prescription of threshold effects at [i = mt, 
while the value m s (Mz) affected by those at /x = m c , /x = mb and n = m t . We consider that 
as for the ratio r 23 at /x = Mz the value is still controversial. Anyhow, we have fixed three 
parameters (ad,oid,€d) only from two values md/m s and m s /mb- 

3.3 CKM mixing 

Next, we discuss CKM quark mixing. Since the parameters (a u ,£ u ) and {ad,ad,£d) have 
been fixed by the observed quark mass rations, the CKM mixing matrix elements |V^ S |, \V c b\, 

|Tub|, an d \Vtd\ are functions of the remaining two parameters <f>i and cj) 2 - In Fig. 2, with taking 
£u = —0.001521, a u = —1.4715, ad = —1.47312, ad = 15.7°, and ^ = 0.004091, we draw allowed 
regions in the (<^>i, fo) parameter plane which are obtained from the observed values [15] of the 
CKM mixing matrix elements and the observed value m of the CP violating phase parameter 
5 q cp in the standard expression of Vckm given by, 

\V US \ = 0.22536 ±0.00061, 

\V cb \ = 0.0414 ±0.0012, 

\V uh \ = 0.00355 ± 0.00015, (3.11) 

l^| = 0.008861“, 

S q CP = 69.4° ±3.4°. 

Here, in order to fix the values of (^>i, fo) we use only two values of the CKM matrix elements as 
input values in the present analysis, so that the remaining tree are our predictions as references. 

As shown in Fig. 1, all the experimental constraints on CKM parameters are satisfied by 


13 


-14.5 

-15.0 

^2 .15.5 

(deg) 

-16.0 

-16.5 

-44 -43 -42 -41 -40 

(j) 1 (deg) 

Figure 1: Allowed region in the (0i, 02 ) parameter plane obtained by the observed values of the 
CKM mixing matrix elements |14 s |, \V cb \, \V ub |, and \Vt d \. We draw allowed regions obtained 
from the observed constraints of the CKM mixing matrix elements shown in Eq. (3.11), with 
taking (a u ,£ u ) = (-1.4715,-0.001521) and (a d ,a d ,£ d ) = (-1.4733,15.694°,+0.004015) . We 
find that the parameter set around (0 1 , 02 ) = (—41.815°,—15.128°) indicated by a star (*) is 
consistent with all the observed values. 

fine tuning the parameters 0 1 and 02 as 

(<£i, 02) = (-41.815°, -15.128°), (3.12) 

which leads to the predicted values for the CKM mixing matrix elements and the CP violating 
phase parameter 5 q CP as follows: 

\V US \ = 0.2261, 

| K cb | = 0.0426, 

\V ub \ = 0.00360, (3.13) 

\V td \ = 0.00920, 

5 q CP = 72.4°. 

In spite of our aim described in the Sec. 1, we are forced to introduce family number- 
dependent parameters ( 0 i, 02 ) hi the present model, too, as the same as in the previous model 
[TO] . However, our aim was to describe all the masses and mixing of quarks and leptons in terms 
of family number-independent parameters expect for the charged lepton masses. Therefore, the 
introduction of the phase parameters 0j were against our aim and unwelcome as it is. In Sec.4, 
we will try to denote these phase parameters 0* in terms of the observed charged lepton masses 
rn ei . 



14 






3.4 PMNS mixing 

Now let us discuss the PMNS lepton mixing. We have already hxed the four parameters 
a u , £ u , 4 >i and as Eqs. (3.7) and (3.12). The remaining free parameters in the neutrino sector 
are only {a vi a v ^ v ). We determine the parameter values of (a„,£j,) as follows: 


(a„, «*„,&,) = (-2.59,-27.3°,-0.0115), 


(3.14) 


which are obtained so as to reproduce the observed values m of the following PMNS mixing 
angles and R u , 


sin 2 2012 = 0.846 ± 0.021, 

, _ = m 2 2 - m 2 ! _ 

Am 2 32 m 2 3 - ml 2 


(7.53 ±0.18) x 10" 5 eV 
(2.44 ±0.06) x 10 -3 eV : 


sin 2 2013 = 0.093 ± 0.008, 

(3.15) 

eV 2 

„ = (3.09 ±0.15) x 10' 2 . 
eV 2 V ' 

(3.16) 



Figure 2: Contour curves of the center, upper, and lower values of the observed PMNS mixing 
parameters sin 2 2012, sin 2 2023, sin 2 2013, and R u in the (a„, a v ) parameter space. We draw the 
curves for the case of £„ = —0.0115 and (0i, <^> 2 ) = (—41.815°,—15.128°) with taking (a u ,£ u ) = 
(-1.4715,-0.001521). We find that the parameter set around {a v ,a, v ) = (-2.59,-27.3°) indi¬ 
cated by a star (*) is consistent with all the observed values. 

In Fig. 2, we show the contour plots of the observed PMNS mixing parameters sin 2 20i2, 
sin 2 2023 , sin 2 2013 , and R u in the (a u , a u ) parameter space for the case of = —0.0115 with 
taking (^ 1 ,^ 2 ) = (—41.815°,—15.128°) and (a u . £ u ) = (-1.4715,-0.001521). It is found that 
all the PMNS mixing parameters are well consistent with the observe values in Eqs. (3.15) and 
(3.16). As shown in Fig. 2, all the experimental constraints on the PMNS mixing parameters 
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are satisfied by fine tuning the parameters a v , a u , and as 

(a v , a„, £„) = (-2.59, -27.3°, -0.0115), (3.17) 

which leads to the predicted values for the PMNS mixing angles, R u , and the Dirac CP violating 
phase parameter 6^ P in the standard expression of Upmns as follows: 

sin 2 2012 = 0.857, 
sin 2 2023 = 0.993, 

sin 2 2013 = 0.0964, (3.18) 

R u = 0.0316, 

S e CP = -76.0°. 

Our model predicts 5 ( CP = —76.0° for the Dirac CP violating phase in the lepton sector. This 
is very interesting because the value shows a size similar to 5 q cp = +72.4° in the CKM mixing 
matrix. 

3.5 Neutrino masses 

We can predict neutrino masses, for the parameters given by (3.7), (3.12), and (3.14), as 
follows 

m v i ~ 0.00046 eV, m v 2 ~ 0.00879 eV, m v3 ~ 0.0502 eV, (3.19) 

by using the input value [15] Am 2 2 — 0.00244 eV 2 . 

We also predict the effective Major ana neutrino mass m (m) in the neutrinoless double 
beta decay as 

(m) = \m u i(U el ) 2 + m u2 (U e 2 ) 2 + m u3 (U e3 ) 2 \ ~ 3.8 x 1CT 3 eV. (3.20) 

In Table 3, we list our predictions of the CKM and the PMNS mixing parameters and quark 
mass ratios and neutrino masses together with the observed values . 

4 VEV relation between P and $ 0 

So far, we have tried to described all Yukawaon VEV matrices (Yf) by using only the 
observed charged lepton masses m e i as input values. We have also tried to understand CP vio¬ 
lating phase only by using phase parameters otf which are phases of family number-independent 
parameters aj. Nevertheless, all such attempts have failed because we always needed a phase 
matrix P in order to fit reasonable CKM mixing and quark mass ratios. In this paper, we accept 
the existence of P, and we try to understand the values of the phase parameters (f>i in P from 
the charged lepton mass values m e j. 

In the present model, we have flavon VEVs with diagonal form, P, P, E, E, 4>o, 4*0; 4> e , 
4> e , and Y e . (Here, we omit “(” and “)” .) In considering combinations of U(3) 8 + 1 scalars 
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Table 3: Predicted values vs. observed values. 



|Ks| 

\V cb \ 

\V ub \ 

\Vtd\ 

$CP 

'12 

'23 

r d 

'12 

r d 

'23 

Pred 

0.2261 

0.0426 

0.00360 

0.00920 

72.4° 

0.0458 

0.0600 

0.0611 

0.0312 

Obs 

0.22536 

0.0414 

0.00355 

0.00886 

69.4° 

0.045 

0.060 

0.053 

0.019 


±0.00061 

± 0.0012 

±0.00015 

+ 0.00033 

- 0.00032 

±3.4° 

+ 0.013 

- 0.010 

±0.005 

+ 0.005 

- 0.003 

+ 0.006 

- 0.006 


sin 2 26*12 

sin 2 26*23 

sin 2 26*13 

Ru [ 10 - 2 ] 

fit 

°CP 

m u i [eV] 

m v 2 [eV] 

m v 3 [eV] 

(m) [eV] 

Pred 

0.857 

0.993 

0.0964 

3.16 

-76.0° 

0.00046 

0.00879 

0.0502 

0.00377 

Obs 

0.846 

0.999 

0.093 

3.09 

- 

- 

- 

- 

< O( 10 “ 1 ) 


± 0.021 

+ 0.001 

- 0.018 

±0.008 

±0.15 







out of those flavons, we have to consider a combination without the parameter A for E and P 
because the R charges of 4>o and <f> e do not contain the parameter A. Only a combination with 
P whose R charge does not include the parameter A is 

(PE ± EP){ = 6j (e^ + ) = 6f 2 cos &, (4.1) 

with R charge of R = |(1 ± A) ± ^(1 — A) = +1. On the other hand, since we have R charges 

i?(<h e ) = 1 . R($ 0 ) = I (4.2) 

for <3? e and 4>o as discussed in Sec.2.3, we have only two combinations which have R charge 
of R = +1, (<E > e )ifc(3 > e) fe ' J and ( < f > o)ia(4 > o) Q ' ? - [Note that (^o^e ± 4> e <l>o) cannot be a candidate, 
because it has R = +1 but it is not a U(3) / singlet.] Therefore, we can take superpotential 

W = Xi[(PE ± EP)Q P \ ± A 2 [(T e <l e + &$ o $o)0p]r (4-3) 

so that we obtain 

k(PE + EP) = $ e $ e + b<!> o^o, (4.4) 

i.e. 

2 k cos <f>i = xj + bxj, (4.5) 

where we have used the dimensionless expressions of P, E, 4>o and 4> e , Eq.(2.9), Eq.(1.7) with 
vo = 1 , and so on. 

Eliminating the coefficient k in Eq.(4.5), we obtain two equations 


COS <fi 1 

xf ± bxf 

(4.6) 

COS (f>3 

1 

± 

cr 

CO to 

COS (f>2 

x\ + bx\ 

(4.7) 

COS (j )3 

x\ + bx\ 
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In Sec.3, we have obtained numerical results fa = —41.815° and fa = —15.128° by putting 
</>3 as (j >3 = 0. In order to avoid confusing, we use notation fa for these numerical results of fa. 
Since we can choose any value of fa in fa —>• fa + fa, we define fa in Eq.(4.5) as 

fa = fa + fa, fa = fa + fa, fa = fa- (4.8) 

The equations (4.6) and (4.7) have two unknown parameters fa and b under the input values 
fa and fa. So, we obtain 

fa = -45.903°, 6 = -1.11586, (4.9) 

which means 

fa = -87.718°, fa = -61.031°, fa = -45.903°. (4.10) 

Regrettably, since we need two input parameters fa and b in order to predict the values 
fa and fa , the present model has no predictability for phase parameters (fa, fa, fa)- (If we use 
the fitting value fa = —41.815° as input value in addition to the input value b = —1.11586, we 
can predict the value fa together with the value of fa.) However, note that the parameters fa 
are family number-dependent parameters, while the parameters fa and b are family number- 
independent parameters. Therefore, the aim of the Yukawaon model that we understand mass 
spectra and mixings of all quarks and leptons only in terms of charged lepton mass spectrum 
and without using any other family number-dependent parameters has been achieved in this 
scenario. 

5 Concluding remarks 

We have tried to describe quark and lepton mass matrices by using only the observed 
values of charged lepton masses ( m e , m T ) as input parameters with family number-dependent 

values. Namely, we have investigated whether we can describe all other observed mass spectra 
(quark and neutrino mass spectra) and the quark- and lepton-mixings (the CKM and the PMNS 
mixings) without using any other family number-dependent parameters. In conclusion, as seen 
in Sec.3, we have obtained reasonable results. Our predicted values are listed in Table 3. 

As seen in Sec.3, we have still used the phase matrix P defined by Eq.(l.lO) in order to fit 
the observed CKM mixing parameters similarly to the past Yukawaon models. However, as seen 
in Sec.4, the most remarkable point of the present paper is that we have succeeded in describing 
the family number-dependent parameters (fa, fa, fa) by the family number-independent pa¬ 
rameters fa and b. Therefore, the main aim in the Yukawaon model have been achieved in the 
present work. However, regrettably, the mechanism proposed has no predictability of the phase 
parameters fa, although it transforms unwelcome family number-dependent parameters into 
family number-independent parameters. The mechanism will be improved in a future version. 

The successful results in the present work suggests the following items: (i) the flavor basis 
in which the charged lepton mass matrix M e is diagonal is more fundamental basis in the flavor 
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physics, (ii) The parameters (m e ,m IJj ,m T ) (i.e. (xi,X 2 ,xz) defined by Eq.(1.7)) are fundamen¬ 
tal parameters in quark and lepton physics. Note that the parameter values are 

extremely hierarchical, while the parameter values ( x\,X 2 ,x 3 ) are mildly hierarchical. Under¬ 
standing of the values of ( x\,X 2 ,x 3 ) will be left to our next task in future. Then, the relation 
(m e + + m T )/( y/m e + + y/m T ) 2 = 2/3 [18] may play an essential role in investigating the 

origin of the parameter values (xi,X 2 ,xs). For reference, we give a trial model on the charged 
lepton mass relation within the framework of the present Yukawaon model in Appendices A and 
B, although this is only a trial one. 

In this model, there are four phase parameters a v , ay and (</>i, </> 2 ). The parameters a v and 
ad play a role in giving mass ratios in the neutrino and down-quark sectors, respectively. The 
parameters which purely contribute to the CKM and PMNS mixing matrices as CP violating 
phase parameters are only (czi>i, </> 2 ) - These parameters can commonly contribute to CKM and 
PMNS mixing matrix, so that those play an essential role in both the predicted values of 5 q cp 
and S^p. It is interesting that, in spite of different values between ad and a v , the results of CP 
violating parameters d^p and &CP take a similar magnitude, 5 q CP ~ ~^cp ~ 70 °. 

In conclusion, it seems to be certain that all of the observed hierarchical structures of 
quark and lepton masses and mixings are commonly originated from the hierarchical values of 
(m e , rrifj,, m T ) which are described by the fundamental parameters ( 21 , £ 2 , 2 : 3 ). Of course, the 
present Yukawaon model has to be still improved with respect to the R charge assignments, 
number of flavons, number of adjustable parameters, CP violating phase parameters, and so on. 
In addition to this, our next task is to investigate the origin of the parameters (xi,X 2 ,xs). 
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Appendix A: Charged Lepton Mass Relation in the Yukawaon Model 

The charged lepton mass relation [18] 


K = 


m e + + m T 



2 

3’ 


(Al) 


is one of the main motives of the Yukawaon model in the earlier stage m- The relation (A.l) 
can be understood from VEV of U(3) 8 + 1 scalar, (<3E> e ) = diag(y/m e , \frnr) as 


Tr[<M e ] 

(Tr[6 e ]) 2 ’ 


(d.2) 


where we have omitted VEV notation “(” and “)” for simplicity. Also, hereafter, for simplicity, 
we denote Tr[A] as [A]. However, in the present scenario of the Yukawaon model, there is no 
8 + 1 scalar <h e , but we have only 6 and 6* scalars <h e and 4+. The purpose of the present paper 
is to understand mass ratios and mixings of quarks and leptons under the given parameters 
(m e , m M , m T ), and it is not to investigate that the origin of the values (m e ,m^,m T ). 

However, in this Appendix, let us try to understand the mass relation (A.l) according to 
an idea suggested in Ref. m • First, let us introduce 8 + 1 scalar <f> e . By using the following 
superpotential: 

W = »[$ e 9 e ] + \ e [($ e E + E$ e )e e ], (A.3) 

we obtain a relation 

<l e = $ e £ + £<f> e . (A.4) 

Since R(E) = |(1 — A) as seen in Eq.(2.19), 4> e has R charge as 

R@e) = 1-^A. (A.5) 

Let us take A = +1, so that we have 

R(4>e) = R($ e ) = (A.6) 


This choice (A.6) causes no problem because 4+ and <f> e have different transformation under 
U(3)xU(3) / . We will comment on the choice R(E) = 0 later. 

Since i?(<h e ) = 1/2, we assume the following superpotential 

W=j (A[<l e <l e ] 2 + A'[4> e ]) 2 [l> 8 4> 8 ]) , (A.7) 


20 




where $8 is an octet part of the nonet <h e defined by 


^8 = ^e-^[$e] I- (-4-8) 

The first term in Eq.(A.7) is the conventional nonet-nonet term. The second term is 
an (octet-octet) x (singlet-singlet) interaction term [19] although the second term is still SU(3) 
invariant. In order to derivate the relation (A.l), the assumption of the second term is essential. 
By noticing that the second term can be expressed as 

[Me][$e] 2 -^e] 4 , (A9) 

we obtain 

= j {2 ( 2 A[d e $ e ] + A'[4> e ] 2 ) $ e + 2A' ^[$ e $e] - ^[4> e ] 2 ) [$ e ]lj • (A10) 

The coefficients of <f> e and 1 must be zero in order to have a nontrivial solution of <f? e (non-zero 
and non-unit matrix form). Thus, we demand 

2A[<M e ] + A'[<i> e ] 2 = 0, (All) 

and 

[$ e 4> e ]-^[4* e ] 2 = 0. (A.12) 

Eq.(A.ll) requires a special relation between A and A'. Note that the relation (A.12) is inde¬ 
pendent of the explicit value of A'. 

Let us comment on the choice of A = +1. This choice means that R(E ) = 0, so that a 
U(3) nonet (EE) takes R(EE) = 0. Therefore, the factor EE can be inserted into any terms 
with R = 2 in the superpotential. However, since (EE) = 1, this does not affect the obtained 
VEV relations practically. The choice A = +1 also gives R charges of Sf as 

(R(S„),R(S d ), R(S e ), R(S U )) = (|, 1, |,2) . (A13) 

It is interesting that the values (1/2,1,3/2, 2) in (A.13) are the same as the values of Ty as seen 
in Eq.(2.30), but the arrangements are different, i.e. (e,u,u,d) for R(Qf), while (u,d,e,u) for 
R(S f ). 
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Appendix B: Alternative Scenario for Charged Lepton Mass Relation 

In Appendix A, we have introduced the new flavon <f> e in addition to the flavons 4> e and 
d» e . So far, we have adhered the idea that the Yukawaon VEV structures take a universal form 
Yf = + £/l (/ = e,u,d,u). However, if we accept an idea that a structure of Y e is 

exceptional, we can introduce ( h e without introducing d» e and <& e as following 

Y e = <M e , <l e = $o^o- (5. 1) 

This expression (B.l) is rather simpler compared with the expression Y e = with P e Q e P e = 

d'oS'e^o given in Sec.2. Therefore, in this scenario, without <f> e and "he [i.e. without Eqs.(A.3) - 
(A.6)], we can use (A.7), so that we can obtain the charged lepton mass relation (A.l). 
However, in this scenario, since we have R charges 

R(Y e ) = + 1, R^e) = \j R(Q o) = \, (B. 2) 

we cannot put the term in Eq.(4.3). In order to avoid this trouble, in a superpotential for 

PE + EP, we a little change the scenario in Sec.4. We assume a mechanism similar to Eqs.(2.1) 
- (2-3): 

W P = [(Ai(PE + EP) + A 2 <Me) Op] + [A [{PE + EP) + A' 2 i e l> e ] [0 P ], (B. 3) 

so that we obtain 

k{PE + EP) = $ e $ e + S P l, (BA) 

i.e. 

2k cos <f>i = xf + £p, (B-5) 

instead of Eqs.(4.3) and(4.4), respectively. From Eqs.(4.5) and (4.6) with instead of b terms, 
we obtain numerical solution 

cj) 0 = 29.222°, = -5.9619, (5.6) 

so that 

0! = -12.623°, 02 = 14.069°, 0 3 = 29.222°. (5.7) 

Maybe, other scenarios are also possible. The purpose in this paper is not to propose 
a scenario which derives the relation (A.l) but to demonstrate a possibility that the family 
number-dependent parameters ( 0 i, 0 2 , 03 ) can, in principle, be described by a family number- 
independent parameter. More reasonable scenario will be given in future. 
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